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Abstract 

Background: Determination of regional lung air volume has several clinical applications. This study investigates the 
use of mid-tidal breathing CT scans to provide regional lung volume data. 

Methods: Low resolution CT scans of the thorax were obtained during tidal breathing in 1 1 healthy control male 
subjects, each on two separate occasions. A 3D map of air volume was derived, and total lung volume calculated. 
The regional distribution of air volume from centre to periphery of the lung was analysed using a radial transform 
and also using one dimensional profiles in three orthogonal directions. 

Results: The total air volumes for the right and left lungs were 1035 +/- 280 ml and 864 +/- 315 ml, respectively 
(mean and SD). The corresponding fractional air volume concentrations (FAVC) were 0.680 +/- 0.044 and 
0.658 +/- 0.062. All differences between the right and left lung were highly significant (p < 0.0001). The coefficients 
of variation of repeated measurement of right and left lung air volumes and FAVC were 6.5% and 6.9% and 2.5% 
and 3.6%, respectively. FAVC correlated significantly with lung space volume (r = 0.78) (p < 0.005). FAVC increased 
from the centre towards the periphery of the lung. Central to peripheral ratios were significantly higher for the 
right (0.100 +/- 0.007 SD) than the left (0.089 +/- 0.013 SD) (p < 0.0001). 

Conclusion: A technique for measuring the distribution of air volume in the lung at mid-tidal breathing is described. 
Mean values and reproducibility are described for healthy male control subjects. Fractional air volume concentration is 
shown to increase with lung size. 
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Background 

X-ray computed tomography (CT) is excellently suited 
to visualisation and accurate measurement of air volume 
in the airway tree [1]. CT scanning provides maps of the 
linear attenuation coefficient of x-rays at an average 
energy of 70 kV. At this energy the attenuation coefficients 
of low atomic number tissues such as soft tissue are closely 
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proportional to density. Thus CT images of the lung reflect 
density and therefore have the advantage of providing 
accurate values of the air volume in the lung. They give 
images of the regional distribution of air volume in 
addition to enabling values of total lung volume to be 
assessed. With multi-detector high resolution CT, the 
whole lung field can be imaged in less than 10 s enabling 
images to be obtained while subjects hold their breath at a 
particular level of respiration. Images at total lung capacity 
are most common, but images at other levels of respir- 
ation can be acquired [2]. Image analysis of the CT 
data allows segmentation of approximately the first six 
generations of the airway tree, enabling accurate definition 
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of the geometry of this section of the airway [3]. This in- 
formation can be applied to derive functional parameters 
using computational fluid dynamics [4]. While the periph- 
eral and alveolated airways are too small to be visualised 
individually, the pattern of their appearance in the images 
can be used to characterise typical patterns occurring in 
different diseases. A simple example is the assessment of 
the density of the peripheral area to give an emphysema 
score [5]. More complex pattern recognition algorithms 
can be used to identify the appearance of different diseases 
[6]. In particular, a method involving parametric response 
maps derived from CT data acquired at full inspiration 
and full expiration has been described, which is useful in 
distinguishing between emphysema and functional small 
airways disease [7]. 

CT-derived parameters have also been used to follow 
developmental and compensatory lung growth [8] and to 
assess structural changes in idiopathic pulmonary fibrosis 
[9] and emphysema [10]. 

Another area where measurements of regional lung 
volume are used is in normalisation of measurements of 
the regional distribution of aerosol deposition using radio- 
nuclide imaging. The relative deposition in conducting 
and peripheral airways is often assessed by a central to 
peripheral ratio. The values obtained depend quite signifi- 
cantly on the size of the areas chosen. This dependency 
can be considerably reduced by normalising the central to 
peripheral ratios to lung volume [11,12]. Lung volume is 
often estimated using transmission scanning or ventilation 
imaging. Both methods have disadvantages - delineation 
of the lungs on transmission is imprecise and ventilation 
does not reflect lung volume in poorly ventilated areas. 
The ability to use CT for this purpose overcomes these 
limitations. The disadvantages of CT are its relatively high 
radiation dose and its accessibility. Dose may be reduced 
using low dose protocols [13] and the use of CT in 
conjunction with both gamma camera tomography 
(Single Photon Emission Computed Tomography, SPECT) 
and positron emission tomography (PET) is becoming 
increasingly commonplace. 

Measurements of lung anatomy from CT have a 
further application in the field of inhaled aerosol de- 
position. They provide anatomical information which 
allows improved computer modelling of the deposition 
pattern. In this respect images of lung volume obtained 
at mid-tidal breathing such as those obtained in some 
SPECT-CT protocols are particularly useful in that 
they provide information on lung anatomy at the same 
stage of breathing as is frequently used in experiments 
on aerosol inhalation. 

An important disadvantage of CT, and indeed a fun- 
damental problem of most three-dimensional imaging 
studies, is that the images are acquired in the horizon- 
tal position. It is known that lung anatomy varies very 



significantly between supine and erect positions; the 
FRC of the lung is greater by around 40% with the 
subject erect compared to supine [13]. Therefore all 
information regarding lung anatomy including lung 
volume values have to be compared in the light of this 
large effect. 

In addition, it has been shown that even when volume 
measurements in the supine position are compared, there 
are significant differences between values from CT and 
other methods, with CT values being generally lower than 
other techniques [14-16]. While CT scanning is now men- 
tioned in standard texts on lung volume assessment [1], 
there is a shortage of reference values for expected lung 
volume values in control subjects. Reproducibility data is 
available for CT lung volumes obtained at total lung cap- 
acity and residual volume [17]. Also, some work has been 
done on reproducibility of lung volume values from CT in 
preclinical studies [13]. 

In this study measurements of lung volume have been 
made using low resolution (and therefore low dose) CT 
during tidal breathing in 11 healthy subjects in the supine 
position. This has enabled mean values of total lung vol- 
ume and patterns of regional distribution to be obtained 
at a mid-tidal breathing position, together with measures 
of inter-subject variability. Each subject was imaged on 
two separate occasions enabling the intra-subject reprodu- 
cibility of the measurements to be assessed. The results 
presented in this paper are derived from data acquired in 
a previously described experiment [18]. The prime pur- 
pose of that work was to provide data for validation of 
computer models of aerosol deposition by comparison 
with experimental measurements using medical imaging. 
The information on lung volume described herein is a 
by-product derived from that image dataset. 

Methods 

Subjects 

The subjects described in this paper were recruited for a 
project to validate computer predictions of aerosol depos- 
ition by comparison with experimental measurements 
using medical imaging [18]. Eleven healthy, never-smoker, 
male subjects, between 20 and 45 years old were studied. 
Subjects had no evidence of respiratory disease and lung 
function tests within the normal range. This included 
being free from the common cold and rhinitis for at 
least four weeks before entry into the study. The study 
was approved by the Local Research Ethics Committee 
of Southampton University Hospital NHS Trust and 
the Administration of Radioactive Substances Advisory 
Committee, and patients gave written consent to partici- 
pate in it. The nature of the study was explained to the 
subjects and each one signed a consent form. Measure- 
ments were made of height, weight, FEV1 and FRC. These 
are shown in Table 1. 
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Table 1 Subject characteristics: height, weight, FEV1 n and 



FRC 


Subject 


Height (cm) 


Weight (kg) 


FEV n (1) 


FRC ml 


POO 


185 


91 


4.82 


3080 


P01 


180 


83 


5.24 


3170 


P02 


174 


56 


4.27 


3370 


P03 


171 


75 


4.68 


3343 


P04 


173 


80 


3.88 


2800 


H01 


174 


82 


4.85 


2360 


H02 


187 


88 


5.57 


4800 


H03 


177 


78 


4.27 


3170 


H04 


186 


87 


5.48 


4170 


H05 


173 


80 


3.76 


3480 


H06 


179 


N/A 


4.98 


2910 



Acquisition of the CT scans 

The CT scans were acquired as part of a combined 
SPECT-CT protocol on a GE Infinia dual head gamma 
camera with the Hawkeye 4 CT attachment (GE Medical 
Systems, Milwaukee, Wisconsin, USA). The participants 
were placed supine on the couch with the whole lung 
field included in the field of view. The CT scan consisted 
of 90 slices with an interslice separation of 4.42 mm. It 
used a voltage of 120 kVp and a current of 1.0 mA. The 
image took approximately 4 mins to acquire and there- 
fore provided a CT image of the thorax at mean tidal 
breathing. The effective dose received by the subjects 
from the CT procedure was 0.8 mSv. The CT images 
have a resolution of the order of 1 mm in the transaxial 
plane, but this is limited in the axial (superior-inferior) 
direction by the inter-slice separation of 4.42 mm. The 
resolution in all dimensions is blurred, as the image is 
obtained as an average over the tidal breathing cycle. 

Segmentation of the lung space 

The CT images were transferred to a workstation running 
the Portable Imaging Computer software (PICS) process- 
ing system [19]. The analysis was carried out using a fully 
three dimensional approach and therefore the first step 
was to convert the data to 4.42 mm cubic voxels. It was 
next necessary to delineate the outer boundary of the lung 
envelope. This is straightforward in principle in that lung 
voxels are generally easily separated from those containing 
tissue by virtue of their different density and therefore 
different Hounsfield number. However in practice the 
algorithms to delineate the lung envelope are not 
straightforward [20,21]. In this study a semi-automatic 
approach has been applied; the method has been de- 
scribed in detail previously [22] and a brief outline is 
presented here. 

The positions of a seed point in the trachea and the 
hila of the right and left lungs were defined manually 



from an interactive visual display of the images. The hilum 
was taken as the point of the first bifurcation of the main 
bronchus, which was readily visible from a coronal view 
for the right lung or a transaxial view for the left. The 
software then proceeded automatically to define regions 
for each lung and the trachea/main bronchi using a 
threshold-based technique [22]. A typical segmentation 
of the right lung from the CT image in a transverse slice 
is illustrated in Figure 1. 

Formation of the lung volume image 

X-ray CT images consist of a map of attenuation coef- 
ficients for photons produced by the x-ray tube. For 
the energy used in CT, the attenuation coefficient for 
the soft tissue voxels in the lung can be assumed to be 
proportional to density. Therefore the density of the voxel, 
Pv, can be calculated using the following equation: 



where I air> I t i SS ue an d I VO x are the image intensity values 
for air, tissue and the voxel under consideration respect- 
ively. The fractional air volume concentration of each voxel, 
(FAVC) can then be calculated as: 

FAVC=(^l-^j (2) 

where p t is the average density of the tissue and blood 
in the lung, which is taken as 1.05 [23]. The air volume 
in the voxel is then FAVC.v where v is the voxel volume. 
It is current practice to use density as the descriptive 
parameter when describing CT scans. However it is 
considered that this new parameter, FAVC, is more 
appropriate to use when considering measurements of 
air volume. The two parameters are simply related with 




Figure 1 Example of the final segmentation of the right lung in 
the region of the right hilum. The outline remains close to the 
lung edge around the lung periphery and includes the vessels 

present in the area of the hilum. 

v> J 



Fleming et al. BMC Medical Imaging 2014, 14:25 
http://www.biomedcentral.com/1471-2342/14/25 



Page 4 of 14 



FAVC being approximately one minus the density value. 
Two dimensional sections through a typical 3D volume 
image derived from CT are shown in Figure 2. 

Calculation of lung volume 

The lung volume is calculated by summing the air 
volume in all the voxels in the lung envelope [14]. This 
simple summation will result in a small underestimate of 
volume due to the partial volume effect. This occurs due 
to the finite resolution of the images which causes voxels 
on the boundary to have increased density and hence 
falsely low estimates of air volume. However this is 
exactly compensated by voxels just outside the lung having 
reduced density and therefore some apparent air content. 
To correct for this effect, the lung envelope was dilated by 
1 voxel (4.42 mm) and the summation of lung air volume 
performed within this dilated volume. Both right and left 
total lung volumes were calculated. 

Determination of mean volume image 

It is often valuable to describe the average appearance of 
a particular organ of the body by forming mean images 
across a range of subjects. Therefore a mean image of 
lung volume was created for the male healthy control 
subjects included in this study. As the lungs in different 
subjects have different shapes and sizes, it was necessary 
to register each image to a standard template. One of 
the subjects with a close to average total lung volume 
was chosen as the template, and each of the others were 
registered to the template using a radial transform based 
on the definition of the lung outline of the subject under 
investigation and the standard template. This was pre- 
ferred to more complicated schemes, such as active shape 
models, for its simplicity. The transform also required 
definition of a single common anatomical location in both 
the subject and template. The hilum of the lung is in many 
respects the natural choice for the centre of a spherical 
transform, as both airways and blood vessels branch out 
approximately radially from this point. However, the 
hilum is either on, or close to, the edge of the lung and 



therefore cannot be used, as nearly half the radial 
paths from it have a very short or zero length within 
the lung. Therefore a point in the centre of the lung O 
(o x , o y , o z ) is chosen, which is related to the hilum 
position H (h x , h y , h z ) such that: 

Ox = (h x + p x ) /2, o y = h y , o z = h z (3) 

where x, y and z represent the left/right, anterior/poster- 
ior and superior/inferior directions respectively. p x is the x 
position of the lateral edge point on the lung at the same y 
and z co-ordinates as the hilum. The transform is illus- 
trated in Figure 3. The purpose is to determine, for each 
voxel in the template shape, the corresponding position in 
the subject shape and then transfer the value of the CT 
image in the subject voxel to the equivalent template 
voxel. Each voxel V t in the template shape is characterised 
by its direction cosines and the fractional radial distance 
(f t = OV t / OE t ) from the central point (O t ) to the 
equivalent extrapolated position on the edge (E t ). The 
line OE s in the subject shape is defined as that passing 
through O s with the same direction cosines as OV t . For 
each voxel along this line the fractional distance to the 
edge (f s = OV s / OE s ) is calculated. The voxel, V s , with 
the value of f s nearest to f t , is found. The CT image 
value in this voxel is transferred to voxel, V t , in the tem- 
plate shape. This process is repeated for all the template 
voxels. Prior to registration the images were expressed 
in terms of the fractional air volume concentration per 
unit space volume, FAVC, and these absolute values 
were preserved in the elastic transform process. 

The process was applied to the right and left lungs 
separately for one of the images obtained from each sub- 
ject and then a combined right and left image produced. 
These 11 combined images were then averaged on a voxel 
by voxel basis to produce the mean image of fractional air 
volume concentrations. As the outer layers of voxels were 
affected by the partial volume effect resulting in a system- 
atic reduction of the values a correction was applied. 
Using successive erosion of the volume outline by three 



( > 



Figure 2 Example 2D transverse, coronal and sagittal sections from a 3D volume image obtained from a CT scan. 
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Subject shape 



Template shape 



Figure 3 Schematic diagram illustrating the radial transform method of elastic registration of the images. The subject's image of 
fractional air volume concentration (FAVC) is transformed to the template shape by finding for each voxel in the template shape V t . the 
corresponding voxel in the subject shape V s . This is based on the direction cosines and the fractional radial distance of V t relative to the centre of 
the transform O t . The voxel in the template shape, V t; is then filled with the FAVC value at V s . 



voxels the average air volume concentration in each was 
voxel layer was found. As expected the outer two layers 
were reduced compared to the third layer. It was then 
assumed that the average FAVC value in all three layers 
was the same and therefore all voxels in each of the outer 
two layers were multiplied by the appropriate factor to 
equalize the average FAVC. 

One dimensional profiles of the average fractional air 
volume concentration were then derived in three orthog- 
onal directions from the mean volume image. 



between left and right lungs and in intra-subject repeat- 
ability assessment, the paired t-test was used. Correlations 
between parameters were performed using linear correl- 
ation. Random differences in intra-subject variability were 
assessed from the standard error of the estimate of linear 
regression between the two sets of measurements, 
expressed as a percentage of the mean. Comparisons 
between measures of variability were performed using 
the F-test. The statistical analyses were performed 
using Microsoft Excel (Redmond, WA). 



Shell analysis of volume 

The radial transform concept was also applied to determine 
the variation of fractional air volume concentration from 
the centre to the periphery of the airway tree. In this case 
the transform was carried out using the hilum as the centre 
of the transform with the lung being divided into ten con- 
centric shells [24]. For each subject the total air volume per 
shell, Vi and mean fractional air volume concentration per 
shell FAVCi were calculated. The mean and standard de- 
viation fractional air volume concentration per shell 
across all subjects were calculated and the intra subject 
reproducibility was assessed. This data was summarized 
as a central to peripheral ratio defined by 

5 

C/P = ^~ (4) 

i=6 

Statistical analysis 

All systematic differences between lung volume parameters 
were assessed using the t-test. For comparing parameters 



Results 

Total lung volume parameters 

The results of measurements of space volume, air volume, 
and fractional air volume concentration are summarized 
in Table 2. The mean values of all parameters were 
significantly greater for the right lung than for the left 
(p < 0.0001). The coefficient of variation of inter-subject 
variability of air volume was significantly greater than for 
the air volume concentration (p<0.01). The coefficients 
of inter-subject variability for all parameters were greater 
for the left lung compared to the right, but in no case were 
the numbers sufficient to give statistical significance. 

The total air volume in both lungs from CT correlated 
significantly with an erect measurement of FRC (r = 0.67, 
p < 0.01, COV 15%). The mean total volume (1900 ml) 
was significantly lower than for FRC (3332 ml). The su- 
pine volume from CT was therefore 43% reduced from 
the erect value obtained by helium dilution. The inter- 
subject coefficient of variation for FRC (19.9%) was 
lower than that for CT volume (31.6%) but this difference 
just failed to achieve significance. Fractional air volume 
concentration correlated significantly with space volume 
(r = 0.78, p < 0.005) (Figure 4). 
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Table 2 Summary of results of measuring global values of lung volume and fractional air volume concentration 





Mean (ml) 


Standard 
deviation (ml) 


Coefficient of 
variation (%) 


Coefficient of variation of 
intra subject reproducibility (%) 


Right lung space volume 


1506 


334 


22.2 


5.3 


Left lung space volume 


1208 


375 


31.0 


5.8 


Right lung air volume 


1035 


280 


27.1 


6.5 


Left lung air volume 


864 


315 


36.5 


6.9 


Fractional air volume concentration right 


0.680 


0.044 


6.5 


2.5 


Fractional air volume concentration left 


0.658 


0.062 


9.4 


3.6 



The coefficients of variation of repeated measurement 
of right and left lung space volumes and air volumes 
were 5.3% and 5.8%, and 6.5% and 6.9%, respectively. 
The correlation of repeat measurements of lung volume 
for both lungs is shown in Figure 5a. The corresponding 
Bland- Altman plot, showing the variation of the difference 
between the measurements and their mean is shown in 
Figure 5c. This suggests that the repeatability in terms of 
absolute volume may increase with lung size but more 
subjects would be required to confirm this. The coefficients 
of variation of intra-subject repeatability for the mean 
fractional air volume concentration were 2.5% and 3.6%, 
respectively. The correlation of repeat measurements is 
shown in Figure 5b and the corresponding Bland-Altman 
plot in Figure 5d. 

Regional lung volume parameters 

Transverse and coronal slices sampled through the mean 
air volume concentration images are shown in Figure 6. 
The full mean volume data is available on the depart- 
mental website. The results of shell analysis of the air 
volume data are shown in Figure 7. This illustrates the 
variation of mean total air volume and mean fractional 
air volume concentration for both the right and left 
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Figure 4 The variation of fraction air volume concentration 


with lung space volume in control subjects for the right (♦) and 


left (■) lungs. 





lungs. The fractional concentration of air in shell 1 assumes 
an intermediate value for both lungs. Since the origin of the 
radial transform is at the hilum, which is taken as the 
centre of the airway where the main bronchus bifurcates, 
some of the voxels of shell one will contain air, while others 
will contain the airway wall surrounding the bronchus. 
These balance out to give an air volume concentration of 
around 0.5. The air volume concentration then decreases 
with shell number due to the presence of major blood ves- 
sels in shells 2 and 3. The fraction of air then increases with 
shell number reaching a fairly constant mean value of about 
0.68 in the outer 4 shells. The concentration is slightly 
higher for the right lung (0.69) than for the left (0.67). 

The mean value of the central to peripheral (c/p) ratio 
of air volume for the left lung of (0.089 +/- 0.013 SD) was 
significantly less than that of the right lung (0.100 +/- 
0.007 SD) (p< 0.0001). The coefficient of variation of 
intra-subject reproducibility for repeated measurements of 
the c/p ratio was 8.5%. 

The results of forming one dimensional profiles of the 
mean volume data set are shown in Figure 8. Figure 8a 
shows the variation of air volume concentration within 
the lung envelope from lateral to medial position. This 
demonstrates the decrease of air concentration medially 
due to the presence of the large blood vessels in the 
medial region. Figure 8b shows the variation of air volume 
concentration within the lung envelope from anterior to 
posterior position. This demonstrates the decrease of air 
concentration posteriorly, presumably due to the effect of 
gravity in the supine position. Figure 8c shows the vari- 
ation of air volume concentration within the lung enve- 
lope from superior to inferior position. There is a general 
decrease from superior to inferior. This may be due to 
compression of the lower parts of the lung in the supine 
position due to pressure from the abdomen contents. The 
more rapid decrease at the lung base could be an artefact 
of the partial volume effect. Although a correction for the 
partial volume effect has been applied this is only on the 
basis of the average effect over the whole outer voxel layer. 
The lung sections at the base become very narrow, and 
are therefore more affected by the partial volume effect 
than the average for the lung. 
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Figure 5 This illustrates the reproducibility of measuring parameters of lung volume. Graph (a) shows the results for total volume and (b) 
for fractional air volume concentration (FAVC) for the right lung (♦) and left lung (■). In each graph, the line is the line of identity. The corresponding 
Bland-Altman plots of difference in measurements against mean of measurements are shown in graphs (c) and (d). 



Discussion 

This study has demonstrated that it is possible to use 
low resolution, and therefore low dose CT, to provide 
quantitative images of lung volume at approximately 
mid-tidal breathing. These depict the fractional lung air 
volume per unit space volume and they enable, among 
other things, the total air volume of the lungs to be 
obtained. A normal range for this has been defined in a 
set of 11 subjects and an assessment of intra-subject 
variability between repeat measurements has been ob- 
tained. A method for producing a mean lung volume map 
has been described and used to provide such a map for 
this group of subjects. This has been analysed to produce 
a description of the regional variation of fractional air vol- 
ume from (i) centre to periphery of the lung using a shell 
analysis and (ii) as air volume profiles in three orthogonal 
directions. The reproducibility of the central to peripheral 
ratio of air volume derived from the shell analysis has also 
been measured. 

The methodology described in this paper is not fully 
automated, and therefore not ideal for routine use. How- 
ever algorithms are now commercially available, e.g. from 
VIDA Diagnostics, USA, which are capable of fully 



automatic segmentation of the lung and determination 
of lung volume. It will be interesting to see whether 
this software is capable of being used with low dose 
images such as those described in this paper and also 
whether the values obtained of lung volume at mid 
tidal breathing are similar. A key aspect of the segmen- 
tation of the lung and the subsequent analysis is the 
definition of the position of hilum of the lung and how 
much of the area around the hilum is included in the 
lung envelope. In this paper we have used the first 
bifurcation of the main bronchus to define the hilum 
and have used a technique involving dilation and erosion 
to fill in the areas in the lung envelope with high density 
including those around the hilum. It would be useful for 
the lung imaging community to make an effort to stand- 
ardise these definitions so that results of different software 
will be comparable. 

The results on total lung volume are now considered. 
As expected, the right lung space volume was found to 
be significantly higher than the left lung, contributing 
53.9% of the total space volume of the lungs. The average 
fractional air volume concentration was also significantly 
higher for the right lung with the result that the fraction 
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Figure 6 This illustrates the mean 3D air volume concentration dataset. Image set (a) shows transverse slices running from apex to base of 
the lung and image set (b) coronal slices running from anterior to posterior. 
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Figure 7 This illustrates the results of shell analysis of lung air volume. The graphs show the variation of (a) the percentage of total lung air 
volume and (b) fractional air volume concentration from the centre of the lung (shell 1) to the periphery (shell 10). Data for the right (♦) and left 
(■) lungs are shown separately. 



of air volume from the right lung is slightly but signifi- 
cantly higher at 54.8%. Consistent with these findings, 
the averaged coefficient of inter-subject variability was 
found to be higher for the air volume than for the 
space volume, 32% compared to 27%. The inter-subject 
variability of fractional air volume concentration was 
considerably lower at 8.6%. 

The average fractional air volume concentration in 
each lung correlated significantly with the total space 
volume of the lung. This is in agreement with the study 
by Gevenois et al. [25], who observed the fractional air 
volume at full inspiration in the supine position correlated 
with total lung capacity in healthy volunteers. This correl- 
ation may be explained by considering the gas exchange 
requirements of the lung. If density were to remain the 
same as body size increases then both the space volume 



and the air volume of the lung would increase proportion- 
ately. Assuming that alveolar volume is the largest compo- 
nent of lung volume, the volume of alveoli will then 
increase in proportion to body size. However, the surface 
area available for gas exchange will only increase propor- 
tional to body size to the power two thirds. Therefore to 
allow for this, the alveolar volume will have to increase at 
a rate which is greater than proportional to body size to 
maintain the required increase in alveolar surface area. It 
should also be noted that the difference in FAVC between 
left and right lungs also followed this pattern with the 
larger right lung also having the higher FAVC. This 
variation will be important to bear in bind when inter- 
preting results of density measurements in the assessment 
of emphysema. Since most diagnostic CT scanning is 
carried out at full inspiration or full exhalation with 
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Figure 8 Variation of fractional air volume concentration within the lung envelope (a) from lateral to medial position, (b) from anterior 
to posterior position, and (c) from superior to inferior position. 
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high resolution CT, it will be important to confirm this 
finding in subjects imaged at these states of inhalation. 

The mean total, (i.e. left plus right), mid-tidal air volume 
in the supine position used for imaging was 1900 ml. This 
is reduced by 43% from the erect FRC value obtained in 
these subjects using helium dilution of 3332 ml. Using the 
standard equation described by Stock and Quanjer [26], 
the mean predicted erect FRC for these subjects was 
found to be 3377 ml. This is very close to the mean value 
found, so this group of subjects can be considered to rep- 
resent a typical adult male population. Ibanez and Raurich 
[27], found the average male supine FRC using He dilution 
to be 2230 ml. Therefore the CT values found in this 
paper are reduced by 330 ml compared to the helium dilu- 
tion technique. One reason why the mean value obtained 
from CT will be lower is that it only describes the air 
volume from the lung hilum onwards. It does not include 
the nasopharynx, oropharynx, trachea and main bronchi. 
All this volume will be included in a helium dilution 
measurement. The mean total volume of this space in 
adult males is 226 ml. The components of this are naso- 
pharynx, 138 ml [26], oropharynx, 61 ml [28], trachea, 
19 ml [29] and main bronchi 8 ml [30]. This can go some 
way to explain the lower value obtained with CT. How- 
ever, this will be balanced by the fact that a CT image 
obtained during tidal breathing might be expected to 
be higher than FRC as it should be equal to FRC plus 
half the tidal volume. Since the mean tidal volume is 
around 500 ml this means that CT values might be 
expected to be around 250 ml higher. This almost 
exactly balances out the underestimation described 
above due to the missing extra-pulmonary volume. 

The above arguments lead us to consider whether 
there is any technical reason for an underestimation in 
the CT measurements. The most likely explanation here 
is due to the partial volume effect resulting from the 
limited spatial resolution of the scans. A first order correc- 
tion for this effect has been applied as described above. 
However this only uses a global correction to the missing 
volume close to the surface of the lung and it may not 
apply accurately to all parts of the lung where local effects 
may be important. For example towards the base of the 
lung, the organ becomes very narrow and in some slices 
there is a thin section of lung volume, where no voxel 
may come above the threshold. The global correction 
applied will allow for this to some extent, but it may not 
be sufficient. This effect will be greater in low resolution 
CT studies such as those described here, particularly as 
mean tidal breathing has been used as compared to a 
breath hold. Investigation of the magnitude of this effect 
using image simulation techniques will be the subject of 
further study. 

The correlation of CT lung air volume with FRC 
was reasonably good (COV 15%) and comparable with 



previously described results [31]. CT imaging might 
be considered as an alternative method for measuring 
FRC, although the large difference in supine compared 
to erect values would need to be taken into consideration 
in interpreting the results. Further studies comparing CT 
and helium dilution should ideally measure a supine FRC 
with helium dilution to enable a direct comparison of 
results. The inter-subject variability of CT lung air volume 
(32%) was not significantly different from that of FRC 
(20%). However it was considerably larger than for FRC. 
Increased number of subjects will enable clarification of 
whether this difference is significant. 

There was good intra-subject reproducibility for total 
lung air volume in these control subjects, with a coeffi- 
cient of variation of repeat measurements on separate 
occasions being 6.5%. This compares favourably with 
previous assessments of reproducibility [13]. The intra- 
subject coefficient of variation of average fractional air 
volume concentration was even lower at 3.1%. However 
this reduced variation might be expected given the lower 
inter-subject variation for this parameter compared to 
total air volume. 

CT scans of the lung give the opportunity of investi- 
gating the three dimensional distribution of air volume 
in the lung at various stages of respiration. Most work to 
date has centred on looking at the number of voxels with 
very high air volume concentration at full inspiration as a 
measure of emphysema [5]. The values obtained for this 
parameter will be dependent on the precise imaging 
protocol [32] and analysis of these low resolution scans of 
the airway tree at average tidal breathing are not likely to 
be useful for this purpose. The elastic registration of each 
image to a standard template has enabled a mean image 
of the distribution of air volume to be obtained. This 
provides a database of the normal distribution of air 
volume against which air volume distributions of subjects 
with suspected lung diseases might be compared using 
techniques such as statistical parametric mapping [33]. 
The method of non-rigid registration described in this 
paper uses a simple global spatial transformation, suitable 
for the relatively low resolution spatial images acquired. 
More sophisticated techniques would be required for high 
resolution images where fiducial markers such as the 
bifurcation points can be identified using commercially 
available software (e.g. Apollo, VIDA Diagnostics, USA). 
These would allow a more precise local deformation to be 
applied. 

Two methods of analysing the spatial distribution of 
air volume in terms of one-dimensional vectors have been 
described in this paper. The first is the radial transform or 
shell analysis technique which allows a description of the 
variation of air volume concentration from centre to per- 
iphery of lungs. The different shells defined in the analysis 
correlate approximately with airway generation number 
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and therefore provide a first order approximation to 
distribution of air volume concentration with airway 
generation. This form of analysis has been used in de- 
scribing the three dimensional spatial variation of inhaled 
aerosol deposition from the centre to the periphery of the 
lung [24]. In both two dimensional [11] and three dimen- 
sional [34] imaging of aerosol deposition, it has been 
found valuable to normalise the shell data to the distribu- 
tion of either ventilation or lung volume. Normalisation to 
ventilation has the disadvantage that it may not be valid in 
diseases of the lung with abnormal ventilation. Normalisa- 
tion to lung air volume seems a more robust method and, 
as CT provides the best imaging of lung air volume distri- 
bution available, it could be regarded as the gold standard 
method for this type of normalisation. The finding that 
the central to peripheral ratio of lung air volume was dif- 
ferent for the left and right lungs shows that normalisation 
of aerosol deposition should be carried out separately for 
each lung. The intra-subject reproducibility for the c/p 
ratio of 8.5% was considered to be reasonably good, 
although more variation was found for this parameter 
compared to measures of the total air volume in the lung. 

The other method of analysis described is the one di- 
mensional analysis of fractional air volume concentration 
in three orthogonal directions. This has provided some 
interesting observations on the trends of this parameter in 
three dimensions in the supine position. The data provides 
a baseline against which other measurements of these pro- 
files can be compared and it will be interesting to observe 
how consistently the variation described in this paper is 
found in future studies. 

There were also differences in the distribution of air 
volume between the left and right lungs. The one dimen- 
sional profiles illustrated the lower values of FAVC for the 
left lung but the variations in the three orthogonal direc- 
tions all showed similar patterns of variation for both 
lungs. The variation of FAVC with shell number again 
showed the consistently lower value for the left lung, but 
the relative difference was much greater for the inner 
shells resulting in a significantly lower central to periph- 
eral ratio for the left lung. This is presumably as a result 
of the relatively higher proportion of major vessels in the 
reduced space volume of the lung envelope on the left 
side, due to the presence of the heart. 

The results described in this paper on lung volume in 
the mid-tidal position are a by-product of a study carried 
out for a different purpose. For this reason not all aspects 
of the experimental design are optimised for lung volume 
measurement. Nevertheless, it is considered that useful 
novel results on the use of CT scanning to measure total 
and regional lung volume have been obtained. The number 
of subjects studied for example was not particularly 
optimised for this study. However from the repeat 
measurements on the 11 subjects studied, reasonably 



precise measures of average parameters of lung volume 
and repeatability have been calculated. In addition statisti- 
cally significant conclusions on differences between pa- 
rameters and correlations have been made. Comparison of 
the mid-tidal volume measurements obtained in the 
supine position with erect measurement of FRC was not 
ideal. In principle FRC would have ideally been carried 
out supine at the same time as the CT scan was being 
acquired, so that direct comparison of the measurements 
could be made. 

Further work required on the spatial analysis of the 
distribution of volume is a description of the fractional 
air volume concentration by lobe and by segment. Since 
it is now possible to delineate lobes directly from CT im- 
ages, and estimate the outline of the segments, lung volume 
parameters can be expressed per lobe and per segment. Ini- 
tial work has shown that airway volume does vary region- 
ally [35] and this will be the subject of future study. 

The determination of regional lung volume at a par- 
ticular point in the respiratory cycle using CT gives the 
possibility of looking at change in lung volume between 
different points of the cycle. The recent development of 
four dimensional CT (4DCT) enables dynamic CT im- 
aging to be carried out over the respiratory cycle [36]. This 
would enable visualization of regional changes in volume, 
leading to high resolution quantitative ventilation imaging. 
This could potentially give important new information on 
lung physiology. 

The potential of magnetic resonance imaging in pro- 
viding information on lung volume needs also to be 
considered, as this provides three dimensional image 
data without incurring any risk due to ionizing radiation. 
Conventional proton MR does not provide images that 
are easily related to air volume, and so are unlikely to 
compete with CT for this purpose. However, measure- 
ments of regional lung volume would be clearly possible 
using hyperpolarised gas imaging [37]. If this were to 
become readily available then it could find an important 
role in this area. 

Conclusions 

In conclusion a method of delineating the lung outline 
from CT images of the thorax has been described. This 
has been applied to low resolution CT scans obtained in 
control subjects at mean supine tidal volume. Average 
measurements of total lung volume and its regional dis- 
tribution are described together with the inter-subject 
variability and the intra-subject variability of repeated 
measurements. The values obtained are reproducible 
and are in good agreement with previous results where 
such a comparison is possible. Fractional air volume con- 
centration is shown to increase with lung size, a result 
which has important implications for the quantitative 
interpretation of measurements of lung density. 
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